3^2 log(3) mylist=c(3,2,4,6) mylist mylist[1] mylist[[1]] x=3+2*I x=3+2*i x x=3+2i x Re(x) Im(x) P=c(.1.2.3.4) P=c(.1,.2,.3,.4) V=seq(1:4) V sum(V*P) Puromycin read.txt("test3702.txt") read.txt("test3702.txt.rtf") read("test3702.txt.rtf") read("test3702.txt") #This is a comment 1+3 4*5 exp(1) runif(10) runif(10) runif(5,6,9) #vector, lists mylist=c(1,5,23,30,3,4,6,4,,) mylist=c(1,5,23,30,3,4,6,4,6,2) sum(mylist) length(mylist) mylist[1] mylist[11] alist=c(F,T,F) alist clist=c("red","orange","pink") clist[2] rep(2,4) seq(-1,1,0.1) length(seq(-1,1,0.1)) mylist ylist=c(1,2,3,4,5,6,7,8) mylist+ylist ylist=c(1,2,3,4,5,6,7,8,9,0) mylist+ylist #dataframes height=c(155,177,164,180) weight=c(65,83,74,115) d=data.frame(height,weight) d d[1,2] d[1,] d[,2] d[1:3,] d[d$height>170,] d$height>170 mylist mean(mylist) sd(mylist) sd? median(mylist) hist(mylist) hist(mylist,col="green") hist(mylist,col="green",xlab="days of month") hist(mylist,col="green",xlab="days of month",ylab="number") hist(mylist,col="green",xlab="days of month",ylab="number",main="MAT 3702") choose(5,2) choose(2,5) #out of 5 choose 2 choose(5,2) # dbinom(x,n,p)=b(x,n,p) dbinom(2,3,.1) pbinom(2,3,.1) x=1; for(i in 1:2){x=x+1} x x=0; for(i in 0:2){x=dbinom(i,3,.1)+x} x # pbinom P(X<=x), dbinom P(X=x) pbinom(2,3,.1,lower.tail=F) pbinom(2,3,.1,lower.tail=FALSE) 1-pbinom(2,3,.1) #P(X<=x), x=310 P(310,527,105) pnorm(310,527,105) #mean=527, sd=105 pnorm(500,527,105) pnorm(500,527,105,lower.tail=F) pnorm(0) #qnorm(p,m,sd)=x st P(X<=x)=p qnorm(.3985,527,105) qnorm(.3985,527,105,lower.tail=F) qnorm(.25) qnorm(.25,lower.tail=F) x=c(0,1) toss=sample(x,10) x toss=sample(x,10,replace="T") toss=sample(x,10,replace=T) toss toss=sample(x,2) toss toss=sample(x,10,replace=T) toss sum(toss)/length(toss) toss=sample(x,100,replace=T) sum(toss)/length(toss) weight=c(.1,.9) coin=sample(x,100,replace=T,weight) coin sum(coin)/length(coin) save.image("~/3702Jan25.RData") x savehistory("~/MAT3702Jan25.Rhistory")